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Abstract 

We consider a linearized inverse problem arising in offshore seismic 
imaging. Following Nolan and Svmes[2"8]. one wishes to determine a sin- 
gular perturbation of a smooth background soundspeed in the Earth from 
measurements made at the surface resulting from various seismic experi- 
ments; the overdetermined data set considered here corresponds to marine 
seismic exploration. In the presence of only fold caustics for the back- 
ground, we identify the geometry of the canonical relation underlying the 
linearized forward scattering operator F, which is a Fourier integral opera- 
tor. We then establish a composition calculus for general FIOs associated 
with similar canonical relations, which we call folded cross caps, sufficient 
for identifying the normal operator F*F. In contrast to the case of a 
single source experiment, treated by Nolan [25] and Felea[5], the resulting 
artifact is i order smoother than the main pseudodifferential part of F* F. 

1 Introduction 

This article deals with a linearized inverse scattering problem considered by 
Nolan and Symes [23]. Acoustic waves are generated at the surface of the earth, 
scatter off heterogeneities in the subsurface and return to the surface. The full 
inverse problem would use the pressure field at the surface to reconstruct an 
image of the subsurface. We instead consider the linearized operator F which 
maps singular perturbations of a smooth background sound speed in the sub- 
surface, assumed known, to perturbations of the resulting pressure field at the 
surface. The goal is to left-invert F; standard techniques suggest studying left 
invertibility of the normal operator N — F*F. To start, we make two as- 
sumptions: (i) no single ray connects a source to a receiver; and (ii) no ray 
originating in the subsurface grazes the surface. Under these assumptions, in 
the case of a single source and receivers ranging over an open subset of the sur- 
face, {X3 = 0}, Rakesh[2Ti] showed that F is a Fourier integral operator (FIO). 
Beylkin lj showed that if caustics do not occur for the background soundspeed, 
F*F is a pseudodifferential operator (^DO). 

*The second author was partially supported by NSF grants DMS-0138167 and DMS- 
0551894. 



For more general data acquisition geometries, the canonical relation of F 
depends on the sets of sources and receivers. Nolan and Symes[2E] proved that, 
if both sources and receivers vary over open and bounded subsets E r and S s of 
the surface, then under the traveltime injectivity condition (TIC), generalizing 
the no-caustic assumption, F*F is still a ^>DO. The same result was stated 
by ten Kroode, Smit and Verdel[18] and their proof was completed by Stolk[30] 
who also relaxed the TIC condition in low dimensions. 

For applications in three spatial variables, an important problem is to un- 
derstand the nature of F and F*F for the marine data acquisition geome- 
try [28], where measurements are made on the codimension one submanifold 
S rjS = {(r%, T2\ Si, s%) £ S r x S s : S2 = r2}. This arises as follows: a seismic 
vessel trails behind it both an acoustic source and recording instruments. The 
point source consists of an airgun which sends acoustic waves through the ocean 
to the subsurface. Reflections occur when the sound waves encounter singulari- 
ties in the material of the subsurface. The reflected rays are received by a linear 
array of hydrophones towed behind the vessel. The vessel then makes repeated 
passes along parallel lines (say, parallel to xi axis). 

The purpose of this paper is to consider the marine geometry under the 
assumption that only the simplest, most prevalent type of caustics, namely 
fold caustics, occur for the background soundspeed. Fold caustics are initially 
defined as follows: A ray departing from a source s in the direction a reaches 
at time t a point denoted x(t, a) in the subsurface. If there is a source s, 
such that the spatial projection map (t, a) — ► x(t, a) has a fold singularity and 
only singularities of this type, then we say that the background soundspeed 
exhibits a fold caustic. By the stability of folds, the maps (t, a) — > x(t, a) 
also have at most fold singularities for all nearby sources s' . However, it seems 
that the natural notion of a fold caustic in the context of the overdetermined 
marine data set considered here is the requirement that the analogous spatial 
projection be a submersion with folds, which is the simplest singularity in the 
non-equidimcnsional setting. This will be elaborated upon in §2 and §4. 

We now introduce the linearized scattering operator F considered in 28J, 18J. 
The model for the scattered waves is given by the wave equation: 

J—^P( x ,t)-Ap(x,t) = S(t)6(x-s) (1) 
p(x,t) = 0, t<0, 

where x G Y = Mj_ = {x G K 3 ,^ > 0} represents the Earth, p(x,t) is the 
pressure field resulting from a pulse at the source s and c{x) is the velocity 
field. The linearization consists in assuming c to be of the form c = cq + 8c, 
where c$ is a smooth known background field. The associated pressure field po 
is also assumed known. The linearization of (fTJ) then becomes 
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1 d 2 5p t . . . 25c(x) d 2 p 

Sp = 0, t < 0, 

where p = po + Now, for a given data acquisition submanifold T, rtS C 
31" x dY and appropriate time interval (0, T), we define the linearized scattering 
operator F : 5c — > sX (o,r)- The assumption (n) ensures that F 1 is an FIO 
( [15] , [18] , [29] , [28] ) and (i) ensures that the composition F*F makes sense. 

In the case of the single source model, with only fold caustics appearing, 
Nolan|25 j showed that F is an FIO associated to a folding canonical relation 
in the sense of [21] (also called a two-sided fold), and stated that the Schwartz 
kernel of the operator F*F belongs to a class of distributions associated to two 
cleanly intersecting Lagrangians in (T*Y \ 0) x (T*Y \ 0). This was fully proved 
in [5]. The corresponding canonical relations are the diagonal A and a folding 
canonical relation, different from the original one, which lies in T*X x T*Y . 

In this article we show that, for the the marine geometry, the linearization 
F is an FIO associated to what we call a folded cross cap canonical relation. 
We then prove that the Schwartz kernel of F*F belongs to a class of distribu- 
tions with a microlocal structure similar to that in the case of the single source 
geometry, but with the order of the non-pseudodifferential part of F*F being 
i lower than in the case of a single source. This means that artifacts arising in 
seismic imaging from the presence of fold caustics are ^ derivative smoother for 
the marine geometry than for the single source geometry. 

Composition of FIOs under other singular geometries arising in integral ge- 
ometry and inverse problems has been previously studied in, e.g., [13] . [9] . [TT] . [12] . 
[25] and [5]. 

The article is organized as follows. In §2 we review some C°° singularity 
theory and define the submersion with folds and cross cap singularities. §3 
is a review of the distribution classes associated to two cleanly intersecting 
Lagrangians, I p,l (Ao, Ai) and the operators which have these as their Schwartz 
kernels. In §4 we show that submersions with folds and cross caps appear 
microlocally in the marine geometry in the presence of the fold caustics, and 
we formulate a general class of canonical relations exhibiting these singularities. 
§5 is dedicated to analyzing a model folded cross cap canonical relation, Co, in 
T*R n x T*K™ -1 ; we establish the composition calculus for F*F, showing that 
F*F E J p ' (A, Co) where Co is a folding canonical relation. Finally, §6 provides 
the extension of this to the general class of folded cross caps. We find a weak 
normal form for any folded cross cap canonical relation C C T*X x T*Y which 
allows us to show that F*F S I p,l (A, C), with C a folding canonical relation in 
T*Y x T*Y. 

We would like to thank Cliff Nolan for the helpful discussions, clarifying 
[25], at the Institute for Mathematics and its Applications, Minneapolis, in 
October, 2005. 
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2 Fourier integral operators and singularity classes 



Let X and Y be manifolds and I m (X, Y;C) denote the class of m-th order 
Fourier integral operators (FIOs), F : £'{Y) — > T>'(X), associated to a canonical 
relation C C (T*X \ 0) X (T*Y \ 0). We will focus on the composition calculus 
for two FIOs. Let C x C (T*X \ 0) x (T*Y \ 0) and C* 2 C (T*Y \ 0) x (T*Z \ 0) 
be two canonical relations and F 1 e r ni (X,Y;d) and F 2 € / m2 (F,Z;C 2 ). If 
Ci x C 2 intersects T*X x At*y xT*Z transversally, then Hormander [IB] proved 
that FioF 2 e I m i+ m i (X, Z; doC 2 ) where Ci o C 2 is the composition of C x and 
C 2 as relations in T*X x T*F and T*Y x T*Z. Duistermaat and Guillemin^ 
and Weinstein[3T] extended this calculus to the case of clean intersection and 
showed that if C\ x C 2 and T*X x Ay.y x T*Z intersect cleanly with excess 
e then A o B e /™i+™2+e/2 ( - Xj z . Ci In each of thege caseSi Ci C2 ig 

again a smooth canonical relation. However, in many interesting problems, these 
assumptions fail, and it is important to analyze the composition and understand 
the resulting operators. It turns out that the geometry of each canonical relation 
and the structure of their projections play an important role. 

Let 7Tl and 7Tr be the projections, to the left and right, from C to T*X and 
T*Y, respectively. If either one is a local diffeomorphism, so is the other one 
and then C, is a local canonical graph. In the case of two canonical relations, 
if at least one of C\ and C 2 is a local canonical graph, then C\ x C 2 intersects 
T*X x At*y x T*Z transversally and the general composition calculus applies. 

Now consider the case when the projections are no longer local diffeomor- 
phisms. When one of the projections is singular, i.e., when the rank of its 
differential is nonmaximal, then the other one is, too, and C is called a singular 
canonical relation. (Note: C is still assumed to be smooth.) 

Although corank(d7Ti)=corank(d7Tfl) at all points, the two projections, itr 
and 7ix, may have similar singularities or quite different ones. The singularities 
considered in this article are folds, submersion with folds and cross caps, which 
we now briefly describe. 

Let / be a smooth function / : V — > W, dim V = dim W = N and 
S := S(f) = {x e V : det(df(x)) = 0}. 

Definition 2.1. / has a (Whitney) fold singularity along S if d(det(df)) ^ on 
S , so that S is a smooth hypersurface, df drops rank by 1 there, and Ker df(x) 
intersects T X S transversally for every i£iS. 

Any map which has a fold singularity can be put into a local normal form: 
f(xi,X2, ■ ■ ■ ,xn) = (xi,x 2 , ■ ■ ■ ,xn-i,x 2 n ) with respect to suitable local coordi- 
nates in the domain and codomain [7]. 

Whitney folds are the singularities denoted by Si t o (in the Thorn theory [7]) 
and by Ei o (in the Boardman-Morin theory [23l [24]) in the equidimensional 
case. The non-equidimensional versions of Whitney folds are submersions with 
folds and cross caps. We note the difference in notation: when dim V < 
dim W, the singularity classes SV.o = £r,o, while, if dim V > dim W, then 
5 r .o = Sr+fc.Oj where k = dim V — dim W. 
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Let / be a smooth function / : V ->• W, dim V = N, dim W = M , N > M. 

Definition 2.2. The map / is a submersion with folds if the only singularities 
of / are of type 51,0, i.e., of type Ejy-M+i,o. 

One checks that / is a submersion with folds as follows. At points where 
rank df > M — 1, by [21] , we can choose suitable adapted local coordinates 
on and W such that / has the form: f(x±, x 2 , ■ ■ ■ , %m-i, %j • • • j£Ejv) = 
(iij^v^M-hAW)- The set S\(f) where / drops rank by 1 is described 
by 5i(/) = {x : §^ = 0, M < i < N}. Then / is a submersion with folds 

if S\(f) is a smooth submanifold, i.e., id(^-\ : M <i< N)\, is linearly 

independent, and if the (N — M + l)-dimensional kernel of df is transversal to 
the tangent space to Si(f) in TV. These conditions can be combined [S3] into 

d 2 h 



det 



dxidxj 



* 0. (3) 

M<i,j<N 



and this is independent of the choice of adapted coordinates. 

There are a finite number of local normal forms for a submersion with folds, 
determined by the signature of the Hessian of / [7] : 



f(xi,X2, ■ ■ ■ ,Xn) = (XI,X2, ■■ ■ , XM-i,X 2 M ± X 2 M+1 ± • • • ± X%). 

In the case relevant here, N = M + 1 and the last entry is a quadratic form in 
two variables, which is either sign definite or indefinite; we refer to these two 
possibilities as elliptic and hyperbolic respectively. 

We now define the second singularity class of interest; like the class of sub- 
mersions with folds, it is stable under small C 2 perturbations. It is now assumed 
that dim V = N, dim W — M with N < M, and g : V -> W is a smooth 
function. 

Definition 2.3. We say that g is a cross cap if the only singularities of g are 
of type Si, , i.e., of type T, lj0 . 

To identify a cross cap, we use the description of [53]. At a point where 
dg has rank > N — 1, we can find suitable adapted coordinates such that 
g(xi,x 2 , ■ ■ .,x N -i,x N ) = (xi,x 2 , ■ • ■ ,x N -i,gi,g 2 , ■ ■ .g q ), where q = M-N+l. 
The set Si (g) where g drops rank by 1 is given by 
Si(g) — {x : =0, 1 < i < q}. Assume that there is an Iq, such that 

d 2 

„ a" (0) 7^ 0. Then, g has a cross cap singularity near if the map \ '■ R W — * R 9 

given by xipx,x 2 , ...x N ) = (§§±, J^-, . . . , satisfies rank d X (0) = q. (No- 
tice that this forces N > q, i.e., M < 2N — 1.) These conditions can be refor- 
mualted as: (i) Si(g) is smooth and of codimension q; (ii) the N x N minors 
of dg generate the ideal of Si(g); and (Hi) Ker (dg) fl TSi(g) = (0). 

As for folds, there is a local normal form for cross caps, due to 
Whitney [52] and Morin[55]: 

g(xi,x 2 , ■ ■ . ,x N ) = (x%,x 2 , . . . ,x N -i,XiX N , . . .x M -nx n ,x 2 n ). (4) 
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3 Distributions and operators associated 
to two cleanly intersecting Lagrangians 

Classes of distributions associated to two cleanly intersecting Lagrangian mani- 
folds were introduced by Melrose and Uhlmann [22 and Guillemin and Uhlmann|14). 
We briefly review their definitions and properties. 

First, one proves that any two pairs of cleanly intersecting Lagrangian sub- 
manifolds are (micro)locally equivalent. Thus, one can consider the model pair 
(A , Ai) where A = T * K ™ = {(x, £) : x = 0} and Ai = N*{x" = 0} = {(x,£) : 
x" = £' = 0} with x' = (xi,X2, • ■ • ,Xf.), and x" = (x k +i,x k +2, ■ ■ ■ ,x n ). One 
defines a class of distributions given by oscillatory integrals whose amplitudes 
are called product-type symbols. Let z — (x, s) be coordinates in R m = R" x R fc 
and (£, <r) the dual coordinates. 

Definition 3.1. S p - l (rn,n,k) is the set of all functions a(z,£,er) g C°°(R" 1 x 
R" x R fc ) such that for every K CC M m and every a € Z" , j3 g Z^ , 7 g Z™ 
there is a c q , ( 3 7 k < 00 such that 



Definition 3.2. [H] Let I p ' l (R n ; A , Ai) be the set of all distributions u such 
that it = ui + U2 with mi g and 



with a € S p '' l '(m, n, k) where p' = p — j + f and Z' = / - |. 

At this point, if X is a manifold of dimension n, we can define the class 
I pl (X; Ao, Ai) for any pair of Lagrangians in T*AT \ cleanly intersecting in 
codimension k. The oscillatory integrals we use are oscillatory integrals in sense 
of H6rmander[14, p. 88]. 

Definition 3.3. [H] u g I p ' l (X; A , Ai) if u = m + U2 + 2 u« where ui g 
I p+l (A \ Ai), m 2 g / p (Ai \ A ), the sum J2 v i is locally finite and v, = Fw l 
where F is a zero order FIO associated to x" 1 where \ '■ T*X \ — > T*R" \ 
is a canonical transformation such that x(Aj) — ^j>J = 0j microlocally, and 



We say that a distribution u £ I r (X; Ao \ Ai) if, microlocally away from Ai, 
u g I r (X; Ao), the standard Hormander class of Fourier integral distributions 
on X associated with Aq. 



Remark 3.4. QI] If u g T p ' l (X; A , Ai) then u g A \ Ai) and also 

ueI p (X;A 1 \A ). 



(5) 




m t g if.'(R n ;Ao,Ai). 



(i 



We will also use the notion of nondegencrate phase functions which 
parametrize two cleanly intersecting Lagrangians, introduced by Mendoza [15]. 
Let A G A n Ai and r C X x R k x (R N \ 0) an open, conic set. 

Definition 3.5. [19] A phase function <j>{x, s, 9) defined on T is a 
parametrization for the pair (Ao, Ai) if 

i) 4>o{ x i 9) := <P(x, 0, 9) where <f>o is a nondegenerate phase function parametriz- 
ing Ao near Ao; and 

ii) (f>i(x, (9, a)) :— cj>(x, , 9) is a nondegenerate phase function parametriz- 
ing Ai near Ao- 

We also refer to </>i(x; 9; a) as a multi-phase function for (Ao, Ai). 

For simplicity, we now focus on the case of codimension 1 intersection rele- 
vant here, i.e., k = 1. Let us consider the following example: If Ao = N*{x' = 0} 
and Ai = N*{x — 0}, with x' = (x 2 , £3, . . . , x n ), then (f>(x,s,6') = x' ■ 9' 
+ X\S I 9' I is a parametrization for (Ao,Ai) since 4>o( x ,9') := 4>(x,0,9') — x'-9', 
which is a parametrization for Ao, and cf>i{x, (9' , a)) := 4>(x, -m, 9) = x' • 9' + xi<7 

is a parametrization for Ai. 



Proposition 3.6. [19] Let pi be a homogenous function of degree 1 such that 
Pi(Ao) = and H P1 (the Hamiltonian vector field associated to pi) is not tangent 
to Ao. If Ai is the flow out from Ao D {pi = 0} by H P1 then there is a 
parametrization cf> for (Ao,Ai) which can be chosen such that ^(x,s,6) — 
Pi(x,d x cf>) and cj>(x,0,9) = </>o with </>o a specified parametrization for Ao. 

Remark 3.7. [15] If (f>(x;9;a) is a multi-function for (Ao,Ai) near Ao G 
Aq n Ai, then, if u E T)'(X) has WF(u) contained in a conic neighborhood of 



where a € S p ' l (X x (R \0)xK), the space of symbol-valued symbols of order 
p, I, defined by: for all K CC X, a £ Z^, £ Z+, 7 G Z%, 

\d^d2a(x;e;a)\<c a0lK (l + \d\ + \a\f-^(l + \a\y-^ (6) 
with p=p+l+E±l - » ; /= -l- 1. 

Finally, we define the classes of generalized (or paired Lagrangian) Fourier in- 
tegral operators, to one of which we will show the normal operator F*F belongs. 
Recall that a canonical relation C C (T*X \ 0) x (T*Y \ 0) is a smooth, conic 
submanifold such that C := {(x, y; £, 77) : (x, £; y, —17) G C} is a Lagrangian sub- 
manifold of [T*(X x y), ujt*(xxY)) 7 i- e j C is a Lagrangian with respect to 

Definition 3.8. If C , C x C (T*X\0) x (T*y \0) are smooth, conic canonical 
relations intersecting cleanly, then I p,1 (Cq,Ci) := I p,l (X, Y; Co, C\) denotes 
the set operators F : £'(Y) — > V(X) whose Schwartz kernels are in I p ' l (X x 



A , then u G T p ' 1 



(X;Ao,Ai) iff 




y ; c 0) c(). 



7 



4 Fold caustics in the marine geometry 



In three spatial dimensions, let s be a fixed source on the surface, 
{.t 3 = 0}; H(x, £) = ^(co(x) -2 — |£| 2 ) the Hamiltonian associated to the smooth 
background soundspeed cq(x) in |(3J); and A s the image of T*R 3 \ under the 
bicharacteristic flow associated to H, which is a Lagrangian submanifold of 
T*R 3 \ 0. The assumption of a (point) fold caustic means that the only sin- 
gularities of the spatial projection tty '■ A s — > Y are folds. We make use of 
the description of A s in Nolan [25]. It can be parametrized by ti nCl the time 
travelled by the incident ray, and the takeoff direction {p\,p 2 ,Pz) € § 2 . We can 
change these coordinates to (sci, X2,Pz) [25] . Hence on A s , x 3 = f(xx, X2,p 3 ) 
and (pi,p 2 ) = (gi{xi, x 2 , P3), 52(^1, %2,Pa))- In this new setting, det dir s = 
■§^{x\,X2,p%) and fold caustics occur where = and 7^ 0. 

In the marine geometry, the source s = (si, S2, 0) is subject to the restriction 
s 2 = f2, so we consider just s\ as an independent coordinate. Fix an s 2 G R, 
i.e., consider a single pass of the vessel, and let A S2 be the union of the flowouts 
{A( Sl . S2 ,o) : si S R}, so that A S2 C T*R 3 \ is an involutive submanifold. We 
say that fold caustics (and no worse) appear for the background sound speed 
if, considering si as a variable, the spatial projection iry ■ A S2 — ► R 3 is a 
submersion with folds. By the structural stability of submersions with folds, this 
condition will then hold for all s' 2 close to s 2 . The presence of fold caustics may 
be characterized as follows. The variables X3 and (pi,p 2 ) are functions of the 
other four: x 3 = f(xi, x 2 , si,p 3 ), (pi,p 2 ) = (gi(xi, x 2 , si,p 3 ), g 2 (xi, x 2 , si,p 3 )). 
The differential diry then becomes: 



dny 

We have 



1 














1 








9f_ 


df 


9/ 


M. 




dx 2 


ds\ 


dp3 



rank diry 



2 if JLL = JL£_ — n 

' 9s 1 9p3 

3, if I^Oorj^O. 



9s 1 / dp3 

Suppressing s 2 , let S± := 5i(7Ty) = {J^ = J^- = 0} be the critical set of 
tty , where rank d7ry drops by 1. At points of S^, Ker (i7ry is two-dimensional 
and spanned by {(0, 0, dsi, Sp3)}. The tangent space to is 

TS A = Ker (d xl , X2>P3 , Sl (^-)) fl Ker <-/,,„,,....; i i. (7) 

where 



5/, , 9 2 / <9 2 / <9 2 / 9 2 / 

and 
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a (^L\-( 92 f 92 f 92 f r . 

a xi , X2 ,p 3 , sl {g Si ) - Iq^Q^ Q X2 Q Sl > Q p3dsi , ds 2>- 

Then, iry is a submersion with folds if 



Si is smooth, i.e., the gradients in ([5]) and §§§ are 
linearly independent, and TSi is transversal to Ker dns, 



if 



a 2 / 
dpi 

d 2 f 



a 2 f 

dsidp3 
d 2 f 



?0. 



(10) 



(11) 



dp 3 ds 1 

Next, we parametrize the canonical relation C of F in terms of si, 2:1,352 
and P3 ; (a%, a 2 , y/\ — M 2 ), the take off direction of the reflected ray, writing 
a = (011,0:2); an d t, the variable dual to time. Following [25j . the canonical 
relation C C r*(£ r , s x (0,T)) x T*R\ is parametrized as 



C = {( si ,ri(-),r 2 (-),t mc (-) +i re /(-), cr (-)>Pi(-)>P2(-)>' r ; 

Zl ,^2,/(-) 1 - T ( c 1 (')"l +.9l(-)) I - T ( c 1 (') a 2 +52('))> 

-r( Co - 1 (.),a) V / T _ H T + P3)} 

where 



/(•) 


= /(aJi,a;2,si,P3); 


ri(-) 


= r i3 -(xi,x 2 ,/(a;i,a;2,si,p3), j = 1,2; 


^mc (' ) 


= ^0(^1)^2)^3); 


W(0 


= t re f(xi,X2,f(xi,x 2 ,si,ps); 


a(-) 


= a(x 1 ,x 2 ,f(x 1 ,x 2 ,s 1 ,p 3 ); 


PiO) 


= p J (x 1 ,x 2 ,f(x 1 ,X2,s 1 ,p 3 ),a),j = 1,2; 


ffiO) 


= gj(xi,x 2 ,si,p 3 ), j = 1,2; and 


co _1 (0 


= c 1 (x 1: x 2 ,f(-)) ■ 



It was proven in [38J that F : £' (R^_) — ► P' (E r .„ x (0,T)) is a Fourier 
integral operator, F S 7?(S r , s x (0,T),R+;C); we now show that the presence 
of caustics of fold type (and no worse) imposes certain conditions on C, namely 
that 7Tr : C — > T*R 3 \ is a submersion with folds and 
7T£ : C — * T*(£ rjS x (0, T)) \ is a cross cap. In fact, with respect with the 
coordinates above, 7r fl : R 7 — > R 6 is given by 



7r fl (a;i,a; 2 ,p3,Si,ai,Q; 2 ,r) = (xi, x 2 , /(■); - t(c 1 (') a l + 5i('))> 

- r(co 1 (-)a 2 +. 92 (-)), 

- r( Co - 1 (.)v / l _ FF + P3)). 
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Thus, 



1 














1 








Of 


df 


df. 


df. 


dx± 


dx 2 




dsi 


At 


A 2 


A 3 


Ai 


By 


B 2 


B 3 


Bi 


Ci 


c 2 


c 3 









-rco" 1 


-tc„ 1 ai 









-tcq 1 

-rc n X «2 



-(c x ai +51) 
-(cq 1 a 2 +32) 



i':sj 








I^OorJ^O 



because the matrix 



-(c/ai + .gi) 
-(c g 1 a 2 +g 2 ) 

is nonsingular [55] • Hence, the critical set 5p := Si(tvr) is a smooth, codimen- 
sion two submanifold. (Recall that by general considerations [4], this must equal 
Si(jtl), and cItyr and ditL must drop rank by the same amount at each point.) 
At these points, Ker dirR = {(0, 0, Sp 3 , Ssi, 5a>i, 5oi2, St)} where Sax, Sa 2 , St 
depend on 5p 3 , Ssi. The tangent space to iSf is 



with 



Ker d XuX 



,df. 

>:2,P3,si,a.i,a.2,T\ q , 



n Ker d Xl . x 



os\ 



,9/ 
dP3 



d 2 f d 2 f d 2 f d 2 f 



dxidpz ' dx 2 dp 3 ' dp\ ' <9si<9p 3 



,0,0,0) (12) 



and 



df 



ix 1 ,x 2 ,P3,si,a 1 ,o 2 ,r( ) — ( 



d 2 .f d 2 f 



dx\ds\ ' dx 2 ds\ ' dp 3 ds\ ' ds\ 



d 2 f d 2 f 

" ^,0,0,0), (13) 



and so we see that Ker d-KR is transversal to TSp because of the condition (fT0| . 
Hence, 7r# is a submersion with folds. Note that without further restrictions on 
/, the projection ttr can either an elliptic or hyperbolic submersion with folds. 

Similarly, with respect to the above coordinates, reordered for ease of display, 
hl : MJ — > R 8 is given by 



ir L (si,xi,x 2 ,ai,a 2 ,p 3 ,T) = (si, n(-), r 2 (-), t inc {-) + t ref (-); cr(-), pi(-), p 2 (-), r) , 
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and thus cIttl = 



/ 


1 
























dri 


dr± 


+ 


dn df 


dri 


+ 


dn df 


dri 


dri 




dsi 


dx\ 




dx3 dxi 


dr-2 


dx 3 dx 2 


dai 


da 2 




dr-2 


Or 2 


+ 


dr 2 Of 


Pro 


+ 


dr 2 df 


dr 2 


dr 2 




ds\ 


dx\ 


dx3 dxi 


dx 2 




dx 3 dx 2 


da i 


da 2 




dt ref 


dt inn 


+ 


dt re f df 


Stint 


+ 


dt ref df 


dt ref 


dt r ef 




ds\ 


dx\ 




dx% dx\ 


dx 2 




dx% dx\ 


dai 


da 2 




da 


da_ 


+ 


da df 


da 


+ 


da df 


da 


da 




dsi 


dxi 


dx^ dxi 


dx2 


dx3 dx2 


dai 


da2 




dpi 


dpi 


+ 


dpi df 


dpi 


+ 


dpi df 


dpi 


dpi 




ds\ 


dx\ 


dx% dx± 
dp2 df 


dx2 


dx 3 dx 2 
dp2 df 


da i 


da 2 




dp-2 


dp2 


+ 


d/>2 


+ 


dp2 


dp2 




ds± 


dx\ 


dx% dx\ 


dx 2 


dx 3 dx 2 


da± 


da 2 


V 


























dl, 





dn df 
dx 3 dp^ 
dr 2 df 
dx 3 dp^ 
. , dtref df 



' dx 3 dp3 
da df 

dx3 dp3 

dpi dpi 
dx3 dp3 
dp2 dp2 
dx 3 dp3 





\ 









1 / 



Since t inc and p% are independent coordinates, we have 



dt, 



dps 



= 0. Also, 



because of the choice of the coordinates x\, X2, %3, it follows that J^- 
at the caustic points [26]. 

f 6, if 



df 

i>:r 2 







From this, it follows that the rank of dn^ 
because the matrix 



7, if 



JLL — JLL = n 

dp3 dsi 
dps i dsi T 1 u , 



/ dri dri dri dri \ 

dxi dx2 dai da 2 

dr2 dr2 dr2 dr2 

dxi dx2 dai da2 

dpi dpi dpi dpi 

dxi dx 2 dai da 2 

. dp2 dp 2 dp 2 dp2 i 

\ dxi dx2 dai da2 / 

is nonsingular [25]. Furthermore, S\{ttl) — S^ = 
and the 7x7 minors of dwL generate the ideal of Sf . Finally, it also follows 
that 



J^j- = 0} is smooth 



Ker dirL = {(0, 5a; i, 5x2, Soil, 5ot2, 5p%, 0)}, 

where 5cti, 5ct2, 5xi, 8x2 depend on 5p%, and thus Ker dw^nTSi 
ttl is a cross cap. 



(0). Hence, 



This leads us to formulate a general class of canonical relations with this 
structure. 

Definition 4.1. Let X and Y be manifolds of dimensions n and n — 1, re- 
spectively, and let C be a canonical relation in (T*X \ 0) x (T*Y \ 0). C is a 
folded cross cap canonical relation if: 

a) ttji : C — > T*Y \ is a submersion with folds, with singular set Si, and 
7Tr(«Si) is a nonradial hypersurface in T*Y \ 0; and 

b) ttl ■ C — > T*X \ has a cross cap singularity along Si and wl(Si) is a 
nonradial submanifold. 

We say that C is an elliptic, respectively, hyperbolic, folded cross cap if ttr 
is an elliptic, respectively, hyperbolic, submersion with folds. 
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Remark 4.2. We will see in §6 (see discussion following Prop. 6.2) that 7T£,(Si) 
is necessarily maximally noninvolutive, defined after (116j) below. 

Here, as usual, a conic submanifold T C T*Y \ is nonradial if for all 
(2A V) £ r, the canonical one-form Vjdyj does not vanish identically on 
XL^T. Since 7Tr(6>i) is the immersed image of Si, ttr(Si) is nonradial at 
ttr(c) iff Yl{ l: *R r li)^* R {dyj) ^ as an element of T*C. We can understand the 
significance of this for the canonical relation arising in the marine geometry as 
follows. Using the fact that df = at the caustics, the expressions for ttr and 
dirR computed above imply that, at c 6 Si, 

3 

^Zi' K R i nj)d^ R {dy J ) = -r ((cq 1 ai + gi)dxi + {c^ ) 1 a 2 + gi)dx 2 ) ■ (14) 

In order to be on T c Si, this must be a linear combination of (fl"2"|) and (JT3J) - 
By dill) , the only possible linear combination is the trivial one; however, by the 
expression for 7Tr, this forces r\ = (0,0,7/3) for some 773 ^ 0. That is, the fold 
caustic surface in Y must be horizontal at this point. While there certainly exist 
background soundspeeds co(-) for which this happens, the most basic examples 
of fold caustics arising from refraction about a low velocity lens (see Nolan and 
Symes [27]) have fold surfaces which are not horizontal. 

Similarly, since ttl\si 1S an immersion, ttl(Si) C T*X is nonradial iff 
YK 7T L,£j)d' !T L(dxj) ^ in T*C. The physical interpretation of a radial point, 
i.e., a point c e Si where this fails, for the marine geometry is less clear and, to 
proceed, we will simply need to assume that such points are absent. 
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5 Model case 



We showed in the previous section that the canonical relation C arising from the 
marine geometry in the presence of fold caustics is a folded cross cap. To help 
understand the nature of the normal operator, we first consider a model folded 
cross cap canonical relation, Co, in (T*M™ \ 0) x (r*M n_1 \ 0), parametrized by 
the phase function 

<hix, y, 9") = (x" - y") ■ 6" + ((x 2 n - ^)y n -i + *„i£_i) *i> 



where {x,y,6") el"x 



pro— 1 



pn-2 



\0), in the region {|#i| > c|0"|}. Here and 



at various points below we use the notation x — (x\, . . . , x n ) = (x", x n -\,x n ) = 
(x 1 ,x"',x n ),y=(y",y n - 1 ) and 6" = {9 u ff"). 

For simplicity, in this section and the next one we will make the choice that 
C is a hyperbolic folded cross cap. This corresponds to the choice of the (— ) 



sign in the (x„ 



x„) term of 4> . There are no significant changes needed in 



the calculations for the elliptic case. 
One easily calculates that 

C = {(x", x n -x,x n ,6" ,-2x n - 1 y n - 1 6 1 ,yl_ 1 6i +2x n y n - 1 6 1 ; 

y", J/n-1,0", - {2x n y n - 1 + (x 2 n - xl_!)) 01 ) 



xi 



: 1^1 >c\6"\, Xi=y u 2<i<n-2, 

- yi + ( x l - x l-i)yn-i + x n y 2 n _ x = o} 



We will verify that the projections of Co to the left and to the right have 
the desired singularities. Note that (x,y n -i,6") are coordinates on Co- Hence, 
reordering the variables for ease of display, 



Tr R (x 1 ,x"',y n - 1 ,e 1 ,e"',x n - 1 ,x n ) = (x 1 + (x^ 



x n— l)Vn— 1 "I" x ny n — \i x > 2/n— 1 

(2x„y„_i + {xl - d i) 



and 







( 1 





A 








D 


D \ 











In-3 



















d-KR - 










1 



























1 





























In-3 














V o 





E 


F 





2x„_i6» 1 


-2y n _ x Q x - 2x n 9 1 j 




where A — 


(x 


2 _ , 
n 




f- 2x 


n Dn — 


uB = 


-2x„_iy 


a -i,D = 2x„y„_i +y 2 n 


-l 


E = -2x n 


h, 


and F = - 


2x n 


Vn-l 


-K 


— x n-l)- 


This is a (2n - 2) x (2n 



matrix, with 

rank ditR 



2n-3, if x„_i = x n + y n -\ = 
2n-2, if x„_i 7^ or x„ + y„_i 7^ 0. 
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Let Si" := Si(ir R ) = {(x,y n -i,0") : x n -\ = x n + y n -i = 0} be the set 
where diTR drops rank by 1. Off of this smooth, codimension two submanifold, 
7Tr is a submersion. The kernel of itr at S^' is spanned by {g| — -, f^— } and 

thus intersects the tangent space T<Sf° transversally. We also note that the 



Hessian, (x n -i,x n 



-(a;„ — xf l _ 1 )8i, is sign-indefinite. We thus conclude 



that ttr is a hyperbolic submersion with folds. 
As for 7Tl, we have 



tt l (x" , x n -t , x n , 61 , 9"' ,y n -i) 
so that 



(a;; 9", -2x n -My n -\9\, (yl_ ± + 2x n y n —i 

)*)> 



( In-2 







V o 





1 





A' 







1 





D' 







1 



B' 
E' 








In-3 












-2x„_i#i 

2(x n + y n -i)0i 



where A' = -2y„_i6»i, B' = -2x n _iy„_i, D' = 2y n _ 1 9 1 , and E' = 2x n y„-i+y*_ 1 . 
This is a 2n x (2n — 1) matrix and 



rank cIttl 



2n-2, if x„_i = x„ + y n _i = 
2n - 1, if x„_i 7^ or x„ + j/ n _i 7^ 0. 



Thus, d7TL drops rank by 1 at <sf°, and off of this set is an immersion. 
The kernel of 7Tl is spanned by — -} and intersects the tangent space 

TS C 

transversally. Also, the rank of the differential of (x,y n -i,9") — > 
(— 2x„_i^i, (2y„_i + 2x n ) 6>i) is 2 at <S>f . We conclude that ttl has a cross 
cap singularity. 

We note for future use the images of Sf" under ir^ and 7Tr: since TTn(Sf") = 
{{a>',-x n ,e",a^e x )} and ir L (S?°) = {{x" , 0, x n , 9", 0, -x^)}, we have 



and 



7TL 



(Sf ) = {x rl _! = = £„ + .X^l = 0} . 



(15) 



(16) 



Note that ir R («Sf °) is a nonradial hypcrsurfacc in T*M™- 1 \0, while n L (S?°) 
is a nonradial, codimension three submanifold of T*M™ \ 0, given by defin- 
ing functions p\ — £ n -i,P2 = £i%n-i,P3 = + with Poisson brackets 
{PI1P2} = l){pi)P3} = 0, {p2,P3} = 0. This is maximally noninvolutive in the 
sense that WT*M"| 7ri ,(Si) nas the maximal possible rank for a codimension three 
submanifold of T*R n \ 0, namely 2n - 4. 
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Next, we calculate the composition CgoCo = {(x,£; y, rf) : 3(z, () such that 
(x, £; z, C) e C$ and (2;, C; y, rf) G Co}. We have (x, £; 2;, C) € C^ iff (z, C; a;, G C 
iff 

Zj = arj, 2 < z < n — 2; 
zi - xi + (z 2 - + ZnX 2 ^ = 0; 

Ci = 6 l<i<n-2; 
Cn = 2z„x n _ 1 £ 1 + x n _i^i; 
Cn-i = -2^n-ia5n-i^i; and 

£n-i = -(4 - - 2z„x„_i£i, 

with (z, £; j/, 77) G Co being determined by the same equations, but where (y, rf) 
replaces (x, £). We thus obtain that Cq o Co consists of (x, £; y, rf) such that, for 
some (z n _i,z„) E R 2 , 

Xi=y%, 2<i<n-2; 
£,i=Vi, 1 < i < n - 2; 
z n -i(x ni - y n -i)£i = 0; 
(x„_i - y„_i)(2z n + x n -i + y n -i) = 0; 
Cn-i = - ((4 - + 2z„x„_i) Ci; 
?7„_i = - ((z 2 - zl_{) + 2z n y n - 1 ) r]i; and 

yi-zi + (z 2 - z 2 _ 1 )(a; Il _i - j/ n _i) + z^x 2 .! - y 2 ^) = 0. 

If x„_i = j/„_i , then xi = yi and = 77^-1, so this contribution to 

Cq o Co is contained in A, the diagonal canonical relation in T*R™ _1 x T*K" _1 . 
If, on the other hand, x„_i ^ y n -i , tben 

Xi=Vi, 2<i<n-2; 
€i=Vh 1 < i < n - 2; 
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l)(3x n _ 1 - Vn-l) t 

Sn-l — 4 ?l! 

+ j/„_i)(3j/„_i - x n _i) 
?7n-i = 1 £i;and 

(aj n -i + y„-i) 2 (a;, l -i - y n -i) n 
xi-yi-\ = 0, 

giving a contribution to Cq o Co contained in Co, where Co is the twisted conor- 
mal bundle, 



ri at* \ i i X n—l + Vn— l) 1 Un~l) n n „ . . . „ 

Go = JV < xi - j/i H = 0, £j - j/i = 0, 2 < i < n - 2 

In conclusion, Cq o Co C A U Co, from which it follows that Co is symmetric, 
i.e., Cq = Co- It is easy to see that both projections from Co have Whitney 
fold singularities. Such canonical relations were introduced in [21] and called 
folding canonical relations; they are also referred to as two-sided folds [8] and 
we will use this latter terminology. One also sees that Co intersects A cleanly 
in codimension 1, and furthermore, A n Co is in fact the fold surface of Co- As 
described in §3, one has a well defined class of distributions associated to the two 
cleanly intersecting Lagrangians (A', C ), namely J p ''(R n_1 x A', Cq), for 

p, I G M. A distribution is in this class iff it has an oscillatory representation, 



/ -iv C^Tl— I - Vn— l)( x n — 1+1/71— l) \^ , f in lil\ rin i I \r 1 

U(X V) = e HOl-2/l 4 ~V )■£ + (a:n-l-l/n-l)Sn-l} 

a{x,y;£";£n-i)d£, 
where a is a symbol-valued symbol, satisfying the estimates: 

with p = p + I + ^,1 = — I — h. 

Next we will show that if F € I m ~i(C Q ) then F*F G I p ' l (A,C ), i.e., 
Kp*p G JP'^R™- 1 x R™" 1 ; A', C ), for some p,iel. We have: 

F/(s) = /" e l{(K ''"^'' ) - e '' +(( ^^"-i ) ^- 1+x "«"~i )ei} a(s,y,0)/(y)^dy, 

where a G S" m+ i Thus, 

^.Hx,„) = /e*<..».«'V,. ( ,^,» o( „, e » )I(Vi ,» ) ^» (( ,». 

After a stationary phase in the variables z", 77", the phase function becomes: 
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<f>(x,y,e",z n -i,z n ) = (x" - y") ■ 9" + (z 2 n - z? l _ 1 )(y n - 1 - x n -i)9i + z n (y 2 i _ 1 - a£_ x )0i 
= (x" - y") ■ 0" + (y„_i - x n -i)[(z 2 n - z 2 n _ x ) + z n (x n -x + y n -i)]6i 

and the amplitude becomes a(x, y, z n -i, z n ; 9") G S 2m+1 . 

Following an idea from [12]. [5]. we now make a singular change of variables, 
T : R" -> I"- 1 , T(0", z„_i,z„) = given by: 

6 = 0<, l<t<n-2, 

£n-l = - {{*n ~ Zn-l) +^n(^n-l + 2/n-l)) 
The kernel of F*F can then be rewritten as 

K F . F (x,y) = [ eW'-*'W'+t*~-™-M"-*h(x,y,Qdt, 
where, using the coarea formula [4, p. 249], 

K x iWi€)= ( a{x,y,z n ^ u z n ;i")-p— 
■/{((^-^-i)9i+*»(»»-i+W»-i)80=-C»-i} 

with J n _i the (n — l)-Jacobian of T and c?^ the arc length measure on 
{((z 2 - z 2 _ ± ) + z n (x n -i + y„-i)) 9i = -£ n _i}. The Jacobian is given by 

Jn-l — |Vz£n-lJ 

= \{2z n -\0\,-(2z n +x n -i +y n -i)9i)\ 

n ( 2 i / x n-l + 2/n— 1 \2 1 l/i | 
= 2 I z „-l + ( z n + ^ ^ J l^l' 

which vanishes to first order at z n _i = z n + = 0. Note that, at 

these points, = £i. Thus, 6 is of order 2m in £" and has a 

conormal singularity of order —1 at — ( x n-i+Vn-i) ^ _ q|. ^ nuS; ^ h as 

an oscillatory representation 

Finally, the kernel of F*F becomes 

K F , F (x,y) = [ e '{^"-^")-C"+(^- 1 -^-i)^- 1 +^- 1 - (3: "- 1+ 4 y "- l)2 h)£}fy X} y . g. p)d frfr 
where one can check that & is a product-type symbol, b £ 5 2 ™ : 1 (2n— 2, n— 1, 1). 
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We now introduce a conic partition of unity in (£, p), with supports in {\p\ < 
2|£|} and {\p\ > resulting in a decomposition Kp*p = K a + K 1 . Letting 



one easily sees that on the region {\p\ < c|£|}, this is a multi-phase function for 
(A',Cq) in the sense of Def. 13. 5i i.e., -0 o (x,y;£) : = i>o(x, y, £, 0) parametrizes 
the diagonal Lagrangian A' and ip(x, y; (£, p)) parametrizes the Lagrangian C' . 
Furthermore, on this region, the amplitude is a symbol- valued symbol, belonging 
to S 2m .-iQR2n-2 x Q^n-i \ ) x M). In view of RemarkEJl K° e F^l"" 1 x 
A', Cg), for some p,l € K. The orders may also be found by applying 
Remark l3~7l to p7| . We see that 

, fn - 1) + 1 2n - 2 1 

p = (2m - 1) + = 2m 

F K ' 2 4 2 

and 

.-(-.)- i-i. 

Hence, this contribution to F*F is in I 2m -h-h (A, <%) := I 2m ~hh (M n_1 ,M n_1 ; A, Co), 
with Co a two-sided fold. 

Next, we show that K 1 e / 2m ~2((7o). First, let s = £ n _i/£i. Then we can 
express K 1 (x 1 y) = J L(x, y, s)ds , where L has the phase function 

y, s; p; £") = (x" ~ y") • C + - Wn-i)ft* + (* - ^zi+^j! ) p 

and amplitude c G S'- 1 - 2m + 1 (R 2Tl - 1 x (R\0) x M"- 1 ). Note that ^(z, ?/, s; p) := 
^(x, y, s; p; 0) parametrizes 

To := 7V*{ S - (a: "~ 1 ^"' l)2 = 0} C T*M 2 "- 1 \ 
and y, s; (p, £")) :=^(x, y, s; (/?,£")) parametrizes 

ri := N* |xi - yi = 0, 2 < i < n - 2, x\ - yi + s(x n - 1 - y n -i) = s - ( x ' n ~ 1 + ^"-^ = g 

Thus, ^ is a multi-phase function in the sense of Def 3.5. Introduce a cutoff 
function, \ G C^(R n ),x = 1 near 0, set 

L°(x,y,s) = J e i *(»*«"''«">x(S)c(a :> i/;«^ , )de"<ift 

andletL 1 = L-£°. WehaveL e I 2m -hi (M 2n - l -T ,T 1 ), L 1 € I 2 " 1- * (R 2n-1 ;ri) 
and WFlL 1 ) is contained in the complement of a neighborhood of Tq. 
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Since K 2 is simply 7r»L, the pushforward of L by tt, the projection tt(x, y, s) — 
(x,y), which is a submersion, we compute the pushforwards of the Lagrangians 
r ,ri. It follows from standard results about pushforwards [Hj that, for u £ 
£' (M 2 " _1 ), 

WF(n*u) Q {(x, y;Z,ri) \ 3{x,y,s,£,r},<T) &WF{u) s.t. 

(x, y) = tt(x, y, s), (|, fj, a) = dw*^, 77)}. 

Using e?7r'(£, 77) = (£, 77, 0), it is then an easy calculation that the push forward 
of To, and indeed any neighborhood of To on which a 7^ 0, is empty, so that 
ir*(L°) e C°°, while 7r*(Ti) C C' Q . In fact, 7r* is a FIO of order -|, and the 
application of 7T* to is covered by the transverse intersection calculus. 

Hence tt* : I r (ri) -> i" r -i(C£), so that 

X 1 = n^L 1 ) mod C°° e /^-^Co), 

and thus K f *f = K a + K 1 e P m ~h>h (A', Cfi. 

In conclusion, we have shown that F*F G I 2m ~2> 3 (A, C Q ). For the single 
source data acquisition geometry, Nolan[5S] and Felea[S] showed that F*F £ 
I 2m '°(A, C), with C a two-sided fold. In that situation the artifact, i.e., the part 
of F*F on C \ A, has, by Remark l3.4l the same strength as on A \ C. However, 
it follows from what we have shown here that for the microlocal model Cq of the 
marine geometry, the artifact is \ order smoother then the pseudodifferential 
operator part. In the next section, we show this in full generality for FIOs 
associated with folded cross caps. 

6 Composition calculus in the general case 

By a well known result of Melrose and Taylor 21j, any two-sided fold can be 
conjugated, via canonical transformations on the left and right, to a normal 
form. Unfortunately, even if folded cross caps could be conjugated to a normal 
form, such a result would presumably be difficult to prove. However, we will 
be able to avoid this problem by finding a weak normal form for the canoni- 
cal relation and thus for a phase function parametrizing it, adapting a method 
originally developed by Greenleaf and Uhlmann [TT1 W2\ for some canonical rela- 
tions arising in integral geometry, for which 7r^ belongs to a class containing the 
submersions with folds and is maximally degenerate. This will be sufficient 
for establishing the composition calculus for general folded cross cap canonical 
relations. 

Theorem 6.1. If F G I m ~i (X, Y; C) is properly supported and C is a folded 
cross cap canonical relation, then F*F E I 2m ~ 2 >i (Y, Y; A, C). Furthermore, 
C is a symmetric, two-sided fold in (T*Y \ 0) x (T*Y \ 0), and AnC equals 
the fold surface in C. 
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Remark 6.2. It was shown in \28f that, for the marine seismic data geometry in 
three spatial dimensions, the linearized scattering map F belongs to 7 1_ 3(S], %S x 
(0,T),R 3 ;C). Thus, in the presence of fold caustics and under the additional 
nonradiality assumptions described below Def. \4-l\ it follows from Theorem \6.1\ 
that the normal operator F*F lies in J2 2 (R 3 ,R 3 ; A, C), with C as above. 

Before establishing a weak normal form for a general folded cross cap, we 
first need to find separate weak normal forms for each of the two projections, 
7Tr and 7T£,. 

Proposition 6.3. Let f : V — ► W be a smooth, conic map, with V a smooth, 
conic manifold of dimension In — 1 and (W,u>w) 0, conic symplectic manifold 
of dimension 2n — 2. Assume that f is a submersion with folds at v$ G V 
and f(V) is nonradial in W . Then, there exist local homogeneous coordinates 
(s",s n -i,a",a"') G R"" 2 xMx (R"- 2 \0)xM 2 onV nearv , and local canonical 
coordinates (y',if) G r*R" _1 \ on W near wq :— f(vo), such that vq — 
(0, 0,e*,0), wo = (0;e|) 7 and 

f(v) = f(s", s^a", a'") = (s", s n ^;a", A s '°(a"' ,a'")) , (18) 

where A' G C°°(IR 2n_:L , S' 2 R 2 *) is homogeneous of degree -1 in ■ and takes values 
in the nonsingular quadratic forms of the same signature as Hess f{vo). 

Proof. In the terminology of |12j . a submersion with folds has clean folds of 
multiplicity one. Prop RTSl is then a special case of [121 Lem. 3.A.1], with slight 
change of notation, and (N, m, n) in [12] being (n — 2, 1, 2) here. □ 

Applying HH) to tt r : C — > T*Y \ and writing a'" = (cr„_i, cj„), we see 
that = Si (ttr) = {<r n -l = <?n = 0} and ir R (S? ) = {i]n-i = 0}. Further- 
more, 

n-l 

w c = ^r(^2 dr]j A dyj) 
i=i 

n-2 

= d(T J A ds J + d (A'^W", o-"')) A ds n _i. (19) 

J'=l 

Hence, 

Ker{uj c \T Sl ) =»• (t^— + ■■■)■ (20) 

Also, as in the case of the model canonical relation, Co, we see that ujc has rank 
2n — 2 on C \ Sf , and has rank 2n — 4 both at <Sp and on <Sp , i.e., restricted 
to TSi . Thus, since uj c = it* L Lu T *x as well, the image itl{Si) C T*A \ 0, 
which is smooth, conic, nonradial and codimension 3, must also be maximally 
noninvolutive. Finally, we note that, as a general fact about canonical relations, 
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for all c = (x ,£o,yo,Vo) & C, the subspace dn L (T Co C) < T {xo ^ ) [T*X) is 
involutive; in particular, for cq G Si, d-Kj,{T Co C) is a codimension 2, involutive 
subspace. We are thus led to establishing a (very) weak normal form for cross 
cap maps into symplectic manifolds reflecting these extra conditions. Note that 
this would apply to a more general class of maps than cross caps, since at this 
point we are only using information concerning the first derivatives. 



Proposition 6.4. Let g : V — ► U be a smooth conic map, with V a smooth, 
conic manifold of dimension In — 1 and (U,u>u) a conic symplectic manifold 
of dimension In. Assume that g has a cross cap singularity, with critical set 
S±. Let vq € Si, no = g{va). Assume that dg(T Vo V) < T Uo (U) is involutive and 
dg (T V0 Si) = T Uo (g(Si)) < T UQ U is maximally noninvolutive for all Vq G S%. 
Then, there exist local homogeneous coordinates 

(t, t) = (t", t n , t", t'") G x R x (R™~ 2 \ 0) x R 2 

on V nearvo, and local canonical coordinates (x,£) G T*R™\0 on U near uq '■= 
g(vo), such that vq — (0,0, e*,0), uq — (0;e*), and, writing t'" = (r n _i,r n ), 

g(t,r) = (t",3„_i(t,T),t n ;r",7„_i(t,T),7„(t,T)) ; (21) 



dg n - 



dr n - 



— 7^ 0; and 



7ik'"=o 



dr 3 



= o = 0, n — 1 < i,j < n. 



(22) 
(23) 



Note that, with respect to these coordinates, Si C {r„_i = r n = 0}, 
E 2 "- 3 := g{S x ) C {x„_i = £„_i = £„ = 0} and 

Ker(cuc\T Sl ) = »• (q^— + ••■)• (24) 



Proof. By assumption, E 2 " 3 C U is nonradial, codimension 3 and maximally 
nonivolutive. By an application of Darboux's Theorem, there exist f |17l Thm. 
21.2.4]) local canonical coordinates (x, £) G T*R"\0 such that S 2 ™" 3 C {x n -i = 
Cn-i = Cn = 0} and u = (0,ej). Letting t 3 = g*{xj), tj = g*(£j) for 
1 5. i __• n — 2, and i„ = <7*(cc„), these functions have linearly independent 
gradients at i>o- If t„-i, t„ are any two defining functions for Si, homogeneous 
of degree 1, then (t, r) := (t", f„, r",r n _i,r n ) are local homogeneous coordinates 
on V, with Si C {r n _i = r„ = 0}. With respect to these coordinates, 

g(t,T) = (t" ) X n - 1 (t ) T),t n ;T",£ n -l(t,T) ) £ n (t,T)) 
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and rank dg{v) = (2n — 3) + rank D(i>), where 



D(«) = 



<9r„ 



(£>„_!,£»„), D n _i,D n eT u U. 



Since g is a cross cap, rank D(u) = l,Vv € Si, so that by interchanging r n _i 
and r n , if necessary, we can assume that D n ~\ 7^ and D„ G K ■ D n -i for 
v € 5i near vo- Furthermore, rotating in the x n -i,£,n-i plane, if necessary, we 
can likewise assume that q^"~* (vq) ^ and ^(v ) = 0, so that 



9r„ 



>> 



<9t„ 



near i;q. 



Let 



by assumption, n(u) < T g t v )U is a codimension 2, involutive subspace. We 
have 

with | ct | >> |6| near vq. A simple calculation shows that a vector 
X = Y^j=i a i~Bx~ + c?f~ e T g ( v -jU is in the symplectic annihilator dg(T v V) u iff 



1 < 3 < « - 2,/3 n = 0, a/3„_i - 6o;„_i - ca,, 



a d^-l ' C 9:ef_i + fl 

gf-} forms a basis for dg(T v V)^. Since 

dg(T v V) is involutive iff dg(T v Vy < dg{T v V), we must have 6(«) = c(w) = 0. 
Thus, ^-\ T '"=o = 0, n— 1 < i, i < n. Relabeling, this finishes the proof of 
PropEll' □ 



Now, let C C [T*X \ 0) x (T*Y \ 0) be a folded cross cap canonical re- 
lation, and apply both PropiO] to tt r : C — ► T*Y \ and PropEl to 
ttl '■ C — > T*X\ near co S C. As in [HI [12], we now attempt to reconcile the 
two resulting coordinate systems on C, (s" , s n -i, o" , &'") and (t" ,t n ,r" ,r"'). 
On TSi (which is the same for both projections), loq has rank 2n — 4. By 
(f20f and ([24]) . since the hypersurface {r ra _i = 0} is transverse to Ker (u>c)> 
it must be locally expressible as {s«-i = F(s",cr")} for some smooth F. 
Now set F(s,er) = -A s > CT (cr'", a"')F(s", a")- then F = tt£/, where /(?/,??) = 
—rj n ^iF(y" 1 rj"). Using (fT9|) . we can find a vector field Y on C which satisfies 

i(Y)cj c =dF = -F(s",a")-^— + 0(\a'"\ 2 ) ■ (da", do"). 

OSn-l 
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Thus, exp(Hf) is a canonical transformation of T*Y \ 0, while exp(Y) is an 
t^c-morphism of C, mapping {r„_i = 0, r" 1 — 0} into {s„_i = 0, a'" = 0}. 
Applying these simultaneously on T*Y and C, respectively, we see that one can 
assume that L := {r n _i = 0, t'" — 0} = {s„_i = 0, a'" — 0}. Restricted to this 
(2n — 4)-dimensional submanifold, u>c is symplectic, so by Darboux there exists 
a canonical transformation $(y", 77") = , <&,,» ) of T*R n ~ 2 \ such that 

*V'k) = *"k and$*(a"| L ) = r"|L. 
Extend $ to # : T*y \ — ► T*Y~ \ by 

^(y" 5 Vn-wn", rj n -i) = i$ y »{y", Ti"),y n -i;<S> v "(y", r?"),??„-i) 
and compose C on the right with the graph of Then, 
7*1 00 = = s",7r|,(a: ri _i),7r|,(a: ri ) = i n ,7i#(y n _i) = s„_i and 7r^(r/') = er" 

form coordinates on C near cq. 

Thus, we have so far shown that if C is a folded cross cap, we may assume 
that [x,y n —i,r]") form (micro)local coordinates on C. Therefore there is a 
generating function S(x ) y n -i ) 6") for C, where S is C°° and homogeneous of 
degree 1 in 8" ([T71 Thm.21.2.18]). Hence, C can be parametrized as 

C = {{x,d x S;do.'S, Vn - 1 ,-6" ',-dy^S)} (25) 

and the phase function x( x i 0") = S(x, y n -U @") ~~ y" ' @" 1S a parametrization 
for the Lagrangian C . We now show that the properties of ttl and 7r# impose 
several conditions on S and its derivatives, forcing the phase function to be very 
similar to the model phase 0o considered in §5. 

To prepare the canonical relation C, we first replicate (fT5|) and (fT6|) . Since 
ttr(Si) is a nonradial hypersurface in T*Y\0, by Darboux's theorem, microlo- 
cally there is a canonical transformation from T*Y \ to T*R" _1 \ taking 
7Tij(<Sf ) to ttr(S 1 °), given by (fT5|) . Similarly, 7T£(5f ) is a codimension three 
submanifold of T*X \ 0, which is maximally noninvolutive in the sense that 
uj t*x\- kl (s^) nas rank 2n — 4. Thus, there exist defining functions Pi,P2,P3 for 
7Ti(5p) with {^1,^2} = 1, {pi,P3} = 0, {p2,P3} = 0. By Darboux's theorem, we 
can find a canonical transformation from T*X \ to T*R ra \ mapping 7T£,(<5>f ) 
to ir L (S?°) given by JT$. 

Using (J25J, 71^(0;, y n _i,0") = {d'eS,y n -i,-0" ',-d yn _^S), so (JXSJ implies 

- rf S/n-i' 5 l{x„_ 1 =0=a„+B„-i} = X n^l ( 26 ) 

From (j2"5)) , ■Ki J {x,y n -\,6") = (x,d x S) and (HHJ) implies 

dx n - 1 S\{ Xn _ 1=0=Xn+yn _jj = (27) 

and 

^x„5'|{a:„_i=0=a;„+i/„_i} = — ^n^i ^1 {x„-i=0=x„+2/„_i } (28) 
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Relation (f27|) means that 



S(x, y n -i,e") = S {x", 9") + xl^Siix, y n -i,6") + {x n + y„-i)S 2 (x, y„- U 9") 
From (j2l)|) we obtain 

<52(£,yn-l>0")l{x„_i=O=x„+j/„_i} = -X n ^!, 

and hence 



(29) 
(30) 



Si{x, y n -i,9") = -x 2 n 0i + x n ^iS 3 (x, y n -i, 9") + (x„ + y n -i)S±(x, y n -i,9"). 



Identity (f3"Tj) implies that 



(31) 



5*2(2;, y n -i, 9")\{ Xn _ 1=0=Xn+ y n _ 1 } — — x n d xl S\[ Xn _ 1 —o— Xn+ y n _ 1 y. (32) 

Now, (|3ll|) and (|52")) imply d xl 5|{ Xra _ 1= o= a ;„+ I/ra _ 1 } = #1, and from (f2"9"| we 
have that d xl S = #1, so S = x x 9\ + S a (x"' ,9"'). At this point, our analysis 
shows that the generating function has the form 



S(x 1 y n _ 1 ,9") = x 1 9 1 + S Q {x"',e"')+xl_ l S 1 [x,y n . l ,6") 

+ (x n + y n -i)[-x 2 n 9i + x n - 1 S 3 (x, y„-i,9") (33) 
+ (x n + y n -i)S±(x,y n -i,9")]. 

Next we consider the differentials dnn and dn^. One computes 



dlTR 



j d x "e"S d Xn _ 1 e"S d Xn e"S 



V d x "y n _ 1 S d Xn _ 1 y n _ 1 S d Xn y n _ 1 S 



dy n _ 1 0"S 
1 



"^J/n-lJ/n-l" 



d e »e"S ^ 


— In-2 



Evaluating at <5>f , by p3|) . this becomes 



/ d x „g„S d e » 9 »S \ 

1 

-J„_ 2 

\ S* 3 -2x„6» 1 +2S* 4 25 4 / 

By assumption, dun has rank 2n — 3 at <Sf . Thus, Sq(x",9") is nonde- 
generate, S , 3|{ a! , l _ 1 =o=a; n +i/ n _i}=0 and S , 4 |{ a; „_ 1= o= :c „+y„_ 1 } = x n 6\. We thus 
have 



S 3 (x,y n -i,0") = x n -xSs(x,yn-i,d") + (x n + y n -x)S & (x,y n -x,6"), (34) 



24 



and 



Si(x,y n -x,9") = x n 9 1 +x n ^ 1 S 7 (x > y n ^ 1 ,6") + (x n +y n _i)Ss(x,y n ^ 1 ,6"). (35) 
Similarly, 



dn 



l\s? 



In- 


2 



































1 



































1 


















dx" x 




d x 


•l 1 ' 


s 






dy r 


-ix' 


s 


de 


"x' 


s 


1 n 

M x"x n 




dx n - 


ix n 




dx n x n . 




dy n - 


-1X„- 


_,s 


d > 


X n - 


_ ± s 


dx" x 


,s 


dx n 


-IX 






s 


d v, 


_ix r 


s 


d t 


"x r 


s 



\ 



Using relation p3|) again, we have 



dir L \ s c 



V 



In -2 

















1 

















1 








x"x" So 











dg" x "So 





2S X 


N 


s 3 








N 


2S 4 - 2x n 9 x 


2Si - 2x n 9 1 






where N = S3 + d Xn _ 1 Si. 

By the cross cap condition, the rank of d7ix|{s n _i=Q=a; n +i/»-i} is 2n — 2. 
Thus, S 3 \{ Xn ^ 1=0=Xn+ y n _ 1 - i = and Si\{ Xn _ 1=0 = Xn +y n _ 1 } = x n 6x, and we see 
that relations (jM)) and (f3"5)) follow from the analysis of ir^ as well as ttr. Putting 
together all the previous relations, the generating function becomes: 



S(x,y n - 1 ,6")=x 1 1 + S {x"',e"') + xl_ 1 S x {x,y n . 1 ,e") 
+ (x n + y n -i)x n y n -i9i 

+ {x n +y n -i)[xl_ 1 S 5 (x,y n ^i,9") (36) 

)S 6 (x,y n - 1 ,9") 
+(x„ +y n -i) 2 S 7 (x,y n -i,9")}. 

Since Sq is nondegenerate, C PI {x n -i = x n = y n -i = 0} is the graph of a 
canonical transformation on T*R n ~ 2 ; we may thus assume that Cn {x n _x = 
x n = y„-i = 0} = {(x",0,9",Q;x",0,9",0)}, so that S {x",9") = x" ■ 9", and 
(1551) becomes: 



S(x,y n -i,0") = x" ■ 9" + x^^Sifayn-uO") + (x n + y n -i)x n y n -i9i 

+ (x n + y n -i)[xl_ 1 S 5 (x,y n -i,9") 

i)Se(x,y n -i,9") 
(x„ + y n -i) 2 S 7 {x. y n -i, 9")] 
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Letting 



M(x,y n ^ 1 ,8") = a;^_ 1 5 5 (x, y n _i, 9") + x n ^ 1 (x n + y n -i)S 6 (x,y n -i,0") 

+ (x n + Vn-l) 2 S 7 (x, Vn-1,0"), 

we now see that the Lagrangian C is parametrized by the phase function 

x(x,y,6") = S(x,y n _ u 9")-y"-0" 

= (x" - y") ■ 6" + {xly^ + x n y 2 n _ x ) 9, (37) 
+a%_ l Si(x,y n -i,0") + (x„ + y n -i)M{x, y n -i, 9"). 

The right projection becomes: 

kr{x" ,y n -i,9" ,x n -i,x n ) = (x" + x 2 n _ 1 dgnSi + (x n +y n _i)dgr>M,y n _ 1 ,9", 

(xl + 2x n y n -i)9 1 + xl^dy^Si + M + (x n + y n -i)d Vn _ x M) 

Since ttr is a submersion with folds, Hess(dnR) : Ker(dn^) — ► Coker(dnR), 
its Hessian at <Sp, is nonsingular. At S x , we have Ker(dTTn) = span{ dx d - , 
and 

HesS ^={ 2sAso 29 1+ 6S 7 \ s o) 

At a fixed point cq € 5f , by a rotation in x n -i, x n , we can diagonalize this 
to obtain 5*6 = at cq. Similarly, 29i + 6S7 represents the d£„ component of 
Hess(jTL); by a canonical transformation, we can assume that this equals 29\ 
at Co, so that SV vanishes at cq. Thus, Sq and <SV are small near Co; note also 
that, by the nondegcneracy of Hessian), d Vn _ 1 S\ + 2S$ ^ near cq. 

We are now ready to prove Theorem 16. II Composing F on the left and right 
with elliptic FIOs of order associated with all of the canonical transformations 
of T*X \ and T*Y \ used above, we can assume that the Schwartz kernel 
of F is represented by an oscillatory integral with the phase function given by 
(|37|) and an amplitude a E S m+ z. 

Let x be the phase function of F*F: 



x = x( z ,y,v") - x(z,x,£") 

I II ll\ II i 2 i 2 2 n / ll\ 

= (Z -y )-T} + Z n yn-lVl+ z nV n - 1 m+ z n-lbl{Z,y n -l,V ) 

+[z n + y n - 1 )M(z,y n - 1 ,r)") - (z" - x")£" - z>„_i£i 

-Z„a£_l£l - Zn-l^l^i^n-liC") - ( z n + X n - 1 )M(z,X n -i,^"). 

We will use stationary phase in z" and rj": set d Z "X — and d^'X = 0, 
where 
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d z "X = V" ~€" + {zn + y n -\)d z "M{z,y n _ l ,r]") - (z n + x n _ 1 )d z nM(z,x n _ 1 ,C) 

+£- 1 (.d x „S 1 (z,y n -i,if)-d x „Si(z,x n - 1 ,Z")) 
d m x = zi-yi+zly n -i+z n yl_ 1 + (z n + y n -i)d r , 1 M(z,y n -i,r]'') 

+zl_ 1 d m S 1 (z,y n - 1 ,r]") 
d m X = Zi-yi + (z n +y n -i)d m M(z 7 y n ^ 1 ,'q")+zl_ 1 d, H S 1 (z,y n ^ l7 ?]") 7 2 < i < n - 2. 

Notice that d\,i v „X is nondegenerate. We may solve these equations implic- 
itly for z" and 77" in terms of the other variables: 

V" = Z" + {z n +x n - 1 )d z nM{z,x n - 1 ,^') - {z n + y n -i)d z nM{z,y n - 1 ,rf') 

-z^_i(0z»Si(z,2/n-i,f7'') - d z » Si{z,x n -±,£,")) 
Zi = Vi - (z n +y n -i)d Vi M(z,y n -i, rj") - z^d^S^z^n-i, rj"), 2<i<n-2; 
zi = yi- zlyn-i - Znvl-i - {z n +y n -i)d Zl M{z,y n - 1 ,r}") - ^ti^Sifz.n,!)"). 

We have that d z »Si(z, y n -i, rj")—d Z "Si(z, x n -i,£") vanishes at y n -i = x n -i 
and 77" = so we can write 

d z »Si{z,y n -i,r)") - d z »Si(z,x n -i,£,") = (j/ n _i - x n ^ 1 )d z//yn _ 1 S 1 (z, y n _ u 77") 

-?')8*W' s ^Vn-irf') 

In a similar way, 

(z n + y n -i)d z "M(z,y n -i,r]") - (z n + x n - 1 )d z »M(z,x n - 1 ,€") 

= {yn-i - x n -i)[d z »M(z,y n -i,ri") + {z n + y n - 1 )d z » yn _ 1 M(z,y n -i,r]")] 
+{v" - £,")(z n + yn-i)d z » v »M(z,y n - 1 ,ri) 

Thus, 77" — £" becomes 

= -(Vn-i — x n -i)[d z "M(z, y n -i, v") 

+ (z n + y n -i)d z "y n _ 1 M(z,y n -i,r)") + zl_ 1 d z » yn _ 1 Si(z,y n -i,r]")] 
x [I + (z n + y n -- L )d z " ri "M(z,y n - 1 ,r]") + zl_ 1 d z > /yn _ 1 S 1 (z, y n _ u t/')] -1 

The phase \ then becomes: 

X = {x" -v")-£" + zl{y n - 1 -x n - 1 )^+z n {y 2 n _ l -xl_ l )^ 1 

HVn-i - a;„_ 1 )[4_ 1 9 !/n _ 1 5 1 (-) + M(-) + (z n + t /n _ 1 )9j /n _ 1 M(-)] 

= - y") ■ C" + (y n -i - x„_i){(z 2 + z n (y n -i +x n -i)) 
+zl- 1 (d Vn _ 1 s 1 (-) + S 5 (-)) 

+{z n + y n -i)[z n -iS 6 {-) + (z n + y n -i)S 7 {-) + 3^ M (•)]} 
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where (•) = (y" , z n -i, z n , y n -i, £") an d the amplitude becomes a e S 2m+1 . 
Let 

N(-) = 4-i(^»-i5i(-) + 5 6 (-)) 

+(^„ + y n _i)[2 n _iS 6 (-) + («« + y n -i)S r (-) + dy^Mi-)} 

and 

Q(-) = Zn^Sei-) + (z n + y n -i)S r (-) + ^„_ 1 M(-) 

with 

dy n -i M = zl_ 1 d Vn _ 1 S 5 + z n - 1 S 6 + Zn-xizn + yn-^dy^^e 

+2(z n + y n -i)S 7 + (z n + y„_i) 2 a WB _ 1 5 7 . (38) 

Repeating the argument from §5, we make a singular change of variables, 
Oi = l<»<n-2, 

On-l = -(46+^(^-l+^-l)?l+iV(')). 

We have 

V Zn _ uZ J n -i = -(2z n _i[9 y „_ 1 5i(-) + 5 5 (-)] + (^ + y„-iR„_ 1 Q(-), 
2z„0! + (z„_i + y„_i)fli + a^iV(-)) , 

so that |V z 6»„_i| = iff 

(z n + y n -i)d Zn _ 1 Q(-) 
* n_1 ~ 2(a„ B _ 1 S 1 (.) + S 5 (-)) 

and 

X n -1 +Vn-1 _ 9 Zn N{-) 

Zn ~ 2 20! ' 

At these points, 



= (^+^-^% 1 --L(^ nA r(.))2 
(5 - Q) 2 

:= {Xn - 1+ 4 Vn - l)2 6 1 + P( 1 /',y n _ 1 ,x n _ 1 ,6"). 
where P(y", y„_i, <c n _i, 9") is 
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£(■)■■= ^(^(O) 2 (39) 

" {Zn + Vn ' l)2 wzjSBt m> ~ {Zn + Vn - l)Q{ ' ] 

pushed forward under the previous change of variables. 

Letting S^°(z, y) denote the critical set in the z, y variables, i.e., for C C 
T*Z X T*Y, note that P(-)\ c a 0/ , = and VP(-)| c c , * = 0, so that 

' v n S l "(z,y) v "S 1 "(z.y) 

^l{x»-i=»»-i} = and VP| {x „_ 1=a „_ l} = 0. (40) 



As in the model case, it follows that KF*F(x,y) has an oscillatory integral 
representation, 

f e i {( ,»-^-0"+(*»- 1 -^-^- 1 +£(e,.- 1 - ( -- 1+ /"- l)2 e,-P)} a{x> V} e> p)dedpj 

where a e S 2m ~ x (2n - 2,n - 1,1). On the region {\p\ < c|£|}, the new 
phase function, ijj(x,y;9; p) is a multi-phase for a pair (A',C") in the sense of 
Def l3.5l ijj(x,y;0;fy parametrizes the diagonal Lagrangian A' and ip(x,y;9; p) 
parametrizes a Lagrangian C . Hence, the contribution to F*F from this re- 
gion is in I p ' l (A, C) for some p,l £ I, and the orders of F*F are computed 
using Remark 3.7 in the same way as for K Q in the model case in §5, so that 
p = 2m — | and I = |. On the other hand, the contribution from {|p| > c|£|} is 
handled in the same way as for K 1 for the model case in §5, giving an element 

Of / 2m -2 (C) C / 2m -3<3(A, C). 

Next, we show that C is a two-sided fold; the fact that C is symmetric just 
follows from A U C = C* o C. We have: 



l A u\ 

y" ,y n -i,6 ,0„_H P + Q-dy n -i p ) ■ 



Xi - Vi + ^-d ei P = 0, 2<i<n-2, 
Z 9 //, (x„_i +y„_i) 2 

- 2/1 ~ ^2 l^n-l 7 ^ " P ) 
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and 



The coordinates on C are (x" , x n -i,y n -i,6i,6"'). Note that p = — (x n -i— y n -i)&i 
L &n-i = t^"- 1 ^"- 1 ) g 1 -|- P; thus, t/j and y\ become: 

y% = Xi - (xn-! - y n -\)d Si P, 2 < i < n - 2, 
(sn-i - y„_i)(a;„_i + y n -i) 2 



yi = xi + 



+ (x„-i - y n -i)d 01 P. 



We now consider the projections 7Tr and ttl- We have 

n R (xi,x"',y n -i, 9i ,9"',x n -i) 

(x n -i - y n -i)(x n -i + y n -i) 2 



- ( 



Xl + 



+ (xn-i +y n -i)d 01 P, 



x'" - (x n -! -y n -i)d0>"P,y n -i;9i,6'" 
(x n -i + y„-i)(3j/ n _i - x„_i) 



-6 



i + P- {x n -! - y n -i)dy n _ 1 P^ 



and 



a;" , x„_i , 6»i , 0"' , y n _i ) 

(x„_i + t/„_i)(3a; n _i - y n -i), 



= yx ,x n -i;V\,v 
One easily computes 

d-lTL 



Ker dirt is spanned by - . 
Similarly, 



fan-i - y n -i)d Xn _ 1 P 



( 4-i 








\ 





1 














1 











J„-3 





V o 






D J 




d x , 




-l — 



P and 



• \ 



/ 1 • • • • 

A B ■ ■ 

1 

1 

I n - 3 

... • -D ) 

where A = -(x n -i - y n -i)d%,„ Xl P and B = J n _ 3 - {x n -i - y n -{]dl,,, x mP. 

Hence, det cItyr — D, as well, and Ker dn R is spanned by g| — -. 



V 
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We have: 

dy^D = -\0i + d 2 yl _P 2d 2 Xn _ iyn _ 1 P + (z n _ x - yn-x)dl n _ iyl _P 

and 

d Xn _ t D = - 2^ n _ Ifc _ 1 P + (ar B _a - Vn-i^^P 

Thus, to show that 7Tr and 717, are both folds, with <5>f = {x n -% = = 
AnC, it suffices to prove that 

1% _ P| + lC- l3 /„-x P l + l^.^l < ^il ( 41 ) 

on a neighborhood of cq; to do this, we need to show that d y2 Pi@x n - X Vn.-x 
and <9 2 2 P satisfy the same estimate near the corresponding point. 

Let T — ^'si(.^fSr(.)) ' From (j39|) . we see that 

= -^ z „7V5 z „ y „_ 1 7V-2(z„ + y n _ 1 )(9j / „_ 1 z„ + l)T 

ZC71 

-{z n + y n _i) 2 9 3/n _ 1 T - (z n + y n -i)d Vri _ l Q - {d Vn _ l z n + 1)Q 

and 

= -^-^„^ nXn _ 1 JV-2(2r n + y n _i)(a ;En _ 1 2; n )T 

Zui 

-{z n + Vn-i) 2 d Xn _ l T - (z n + y n -i)d Xn _ l Q - {d Xn _ 1 z n )Q. 

Thus 

-2(d yn _ lZn + 1) 2 T - 2(z„ + y„_i)(^_ iZn )T 

-4(Z„ + yn-O^n-x^n + 1)^^ - (z n + y^xfd^ T 

-2(d Vn _ lZn + l)d Vn _ x Q - (z n + 2/n-i)^2_ x Q - (^2_ i z„)Q. 

Similarly, 

-2(a,„_ 1 z„) 2 T - 2{z n + yn^id^z^T 

-4(z„ + y n -i){dx n -iZn)d Xn -iT - (z n + y n _i) 2 d 2 z_T 

-2{d Xn _ 1 z n )d Xn _ l Q - {zn + y n -\)d%a_^Q- z n )Q 
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and 

d 2 P = — — (d 2 N)d 2 N —9, Nd 3 N 

"l/„_ll„_l 29\ ^ ZnVn — 1 > ZnXn-l 26\ " Z " Vn - 1 %n - 1 

-2(a x „_ 1 z n )(a a „_ 1 z n + 1)T - 2(z n + y n _ 1 )(6g n _ lWn _ 1 z Il )T 

-2(z„ + y n -i)(d Vn _ 1 z n + l)d Xn _ 1 T - 2(z n + y n -i){d Xn _ l z n )d yn _ l T 

~{z„ + y n -i) 2 d 2 n iXn _T - (d Xn _ 1 z n )d Vn _ 1 Q 

-(z n + y n -i)dl n _ lXn _ x Q - {^ n _ XXn _ x z n )Q - (d yn _ 1 z n + tfdx^Q. 

We have 

d Zn N = z 2 n _ 1 (d 2 Zn y n . 1 S 1 + d ^S 5 ) 

x)[z n -xd Zn SQ + S 7 + ( )d Zn S 7 + d 2 n _ lZn M] + Q. 

from which it follows that 

d Ly n -i N = 2 ^n-l^„_ 1 2„-i(^ n!/Ti _ l 5 1 + 8 Zn S 5 ) + zl^dy^id^^Sx + d Zn S 5 ) 
+ {dy n -i z n + l)[z n -ld Zn Se + S 7 + ( )d Zn S 7 + d 2 n _ iZn M] 

+(z n + yn-i)dy n _ 1 [z n -id z „S 6 + S 7 + (z n + y n -i)d Zn S 7 + d 2 n _ iZn M] 
+dy n -xQ (42) 

and 

d 2 ZnXn _,N = 2z n - 1 d Xn _ 1 z n - 1 {% nVn _ 1 S 1 + 8 Zn S 5 ) + z'-i^CC^-A + d ^S 5 ) 
+ (d Xn _ 1 z n )[z n _ 1 d Zri S 6 + S 7 + ( )d Zn s 7 + d 2 n _ 1Zn M] 

+(z n + y n -i)d Xn _ l [z„-xd Zn S 6 + S- 

+(z n + Vn-i)d Zn S 7 + d 2 yniZn M] + d^Q. (43) 

Similarly, 

d Vn -iQ = d Vn _ 1 z n - 1 S 6 + z n ^ 1 d Vn _ 1 S 6 + {dy n _ 1 z n + l)S 7 

+(z n + y n -l)dy n -,S 7 + d 2 yl _M (44) 

and 

d Xn _ 1 Q — d Xn _ 1 z n -iSQ + z n -id Xrl _ 1 SQ 

+ (d Xn _ 1 z n )S 7 + (z n + y n -i)d Xn _ 1 S 7 + d 2 nlXnl M. (45) 



Using (f38|) . we have 



d y 2 _ x M = ^ z n-i{dy n _ 1 z n ^ l )dy n _ 1 S 5 + 2„_i9 J/ 2_ i S' 5 + d Vn _ 1 z n - 1 S 6 
+z n -idy n _ 1 S 6 + d Vn _ 1 z n - 1 (z n + y n -i)dy n _ 1 S 6 
+z n -i(dy n _ 1 z n + l)d Vn _ 1 S 6 + z n -i(z n + y n -i)d 2 2_^S 6 
+2(d Vn _ 1 z n + 1)S 7 + 2(z n + y n - 1 )d Vn _ 1 S 7 (46) 
+2(z n + y n -i)(d Vn _ lZn + l)d Vn _ x S 7 + (z n + y n ~i) 2 d 2 2 S 7 
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and 

d y n -ix n -! M = 2z„_i(5 Xri _ 1 z„_i)9 a7i _ 1 S , 5 + zl_ 1 dy n _ lXril S 5 + d Xn _ 1 z n ^iS e 

+ z n-l(d X n-i z n)dy rl -iS() + z n-l( z n + 2/n-l)5y„_ ia ;„_i ^6 

+2(a,„_ 1 z„)5 7 + 2(z„ + yn-^d^Sr (47) 
+2(z„ + y„-i)(9 a;7l _ 1 2; n )9 an _ 1 S'7 + (z„ + y n -i) 2 dy n _ lXn _ 1 S 7 . 

Similarly, 

dyn-iZn^ = Z n-l^n-i2»^5 + z n-\9 Zn SQ + Z n _idy nl SQ 

+z n -x(z n + y n -x)dy n _ lZn S 6 + 2SV + 2(z„ + y n -i)d Zn S 7 
+2( ) 2 d 2 yn _ 1Zn S 7 . (48) 

Using - igSJ) and the fact that d Zn N,Q,d Vn _ 1 M, S 6 and S 7 are all 
small near Co, we obtain that 9y n _ ian _ 1 M, <9^ 2 M, c^^M, d Xn _ 1 Q, 

dy n -iQ, d Ly n -i N and d l n x n -i N are as welL Thus » 9 9 2 ^ and 

d Xn _ lVn l P are also small, from which it follows that (I41|) holds, and thus C is 
a symmetric, two-sided fold, with A n C = <Sp . 

Remark 6.5. An interesting question, particularly relevant for the marine seis- 
mic imaging problem, is whether (under an ellipticity assumption on F), F*F 
can be inverted, at least modulo operators mapping H s — > H s ~ 2m + S f or some 
6 > 0. Note that, by Remark 13^41 the order of F*F on C \ A is | less than 
on A \ C . However, the standard technique of parabolic decomposition (see 
[10] ) only results in a decomposition F*F = T x + T 2 with T x e /?™(A) and 

2 ' 2 

T 2 € lf"\ (C) . Since C is a two-sided fold, there is further loss of I derivatives 
in terms of Sobolev mapping properties [21] . We will return to this question in 
i- 
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